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Abstract 

For any spherically symmetric black hole spacetime with an ideal fluid source, 
we establish a dual fluid system on a hypersurface near the black hole horizon. 
The dual fluid is incompressible and obeys Navier-Stokes equation subject to 
some external force. The force term in the fluid equation consists in two parts, 
one comes from the curvature of the hypersurface, the other comes from the 
stress-energy of the bulk fluid. 



1 Introduction 

Holography is perhaps the most important keyword in the area of high energy theory 



CO 
> 

in 

■^j- ; over the last twenty years. In spired by the area law of black hole entropy, 't Hooft 

[1] and Susskind [2] respectively proposed the so-called holographic principle in the 

early nineties of the last century. Since then, this concept is greatly pushed forward 

^ ! due to Maldacena's discovery of AdS/CFT duality [3]. Holographic dual of both string 

theory and gravity has been studied extensively, leading to new developments such 

as AdS/CMT, AdS/QCD, etc., and applications of this idea in areas like holographic 

phase transitions, etc., have now attracted serious attentions from the condensed matter 

H ■ community. 
c3 ■ 

Generically speaking, holographic duality relates gravity in the bulk to certain sys- 
tems with other type of interactions (e.g. gauge interactions) on the boundary. In 
some cases the microscopic detail of the boundary degrees of freedom may be difficult 
to determine or is unimportant/beyond the scope of concern. In such circumstances, 
the boundary degrees of freedom may be coarse grained and analogized by some fluid 
component. This is why gravity-fluid duality comes into play. The first relationship 
between gravity and fluid was found quite long ago [4]. In the subsequent 40+ years, 
similar idea pops up from time to time until very recently gravity-fluid duality becomes 
popular and is considered as a new type of holographic duality [5]- [19]. In some sense, 
gravity-fluid dual may be thought of as a kind of "poor man's holography" [12] because 
of lack of information on the boundary. On the other hand, it also helps to understand 
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the holographic dual of gravity in general, i.e. not to be limited to gravity theories 
with AdS asymptotics. 

Although quite a lot of research papers on gravity-fluid duality have been published, 
most of the works are on a case by case basis. It is shown in [6] that given a particular 
boost /scaling for Rindler spacetime, there exist a conserved stress tensor on boundary 
which imply the Navier-Stokes equation of the dual fluid, more relevant works can 
be found [7]- [10], even for higher curvature gravity case we can also get the same 
conclusion [11] [13]. Another way of realizing gravity-fluid duality is to take some 
special kind of boundary conditions to reduce the dynamical degrees of freedom on the 
boundary of some prescribed spacetime background [14] [15]. To understand gravity- 
fluid duality more thoroughly, it seems necessary to realize gravity-fluid duality on 
more general footing, for instance, removing either the metric dependence or the use of 
specific boost / scaling technique would greatly increase the usefulness of such duality 
relations. Moreover, it is desirable to consider the holographic dual of gravity with 
matter source, because it has long been known that Einstein equation with matter 
source is equivalent to Navier-Stokes equation for the matter source regarded as a fluid 
[16]. If one can establish holographic dual of sourced gravity in general, then it is a 
simple step to establish the duality between the fluid equations of the matter source 
and that of the boundary dual. Since this holographic relationship between the bulk 
and boundary fluid is established via gravity-fluid duality, it may be fair to call it as a 
gravity-mediated holography for the bulk and boundary fluids. The present work is a 
first attempt in establishing this kind of fluid-fluid holography. 



2 Setting up the problem 

Our ultimate aim is to establish a duality relationship between the bulk fluid and the 
boundary fluid. In order to make this idea work, it is necessary to consider a gravity 
theory with matter source. We shall deliberately keep the system as general as possible, 
though there are still some places where certain kind of assumptions/restrictions are 
still needed. 

We begin in an (n-|-2)-dimensional spacetime with coordinates = (t, r, 9 l , ■ ■ ■ , 9 n ). 
The Einstein equation with an an ideal fluid source reads 

Gfj, u = R^u — -Rg^ u = —Ag^ u + T^ u , (1) 

= (p + p)u^Uu + P9^, (2) 

where A is the cosmo logical constant, p, p are respectively the energy density and 
pressure in the bulk, is the proper velocity, and we have set Einstein's gravitational 
constant K — 1. 

We assume that there exists a spherically symmetric black hole solution for eq.(l). 
In Eddington-Fenkelstein coordinates, the line element takes the form 

ds 2 = -f{r)dt 2 + 2dtdr + rW, (3) 
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where f(r) is a function of r only, dQ 2 is the line element on a unit n-sphere which can 
be written as 



dtt 2 = dOl + sin 2 °i(d6 2 2 + sin 2 2 (d0^ + •••))• ( 4 ) 

We do not need to solve eq.(l) to get the solution, just the existence of such a solution 
is sufficient to facilitate our discussions. 

Now consider the timelike hypersurface S c determined by taking r = r c constant. 
The coordinates on the hypersurface can be naturally taken to be x a = (t, 9 1 , ■ ■ • , 6 n ). 
The embedding of this hypersurface in the spacetime equips it with a metric (the first 
fundamental form) h^ iu , which is related to the spacetime metric via 

hfj,p — (jfj,u Tifj,n v , (5) 

where is the unit normal vector of E c . For the spacetime metric given by the line 
element (3), we have 

n,= fo,l=,0,... ,o) , n*=C^,y/f,0,---,0j, n 2 = 1, (6) 

and the line element corresponding to reads 

ds 2 n+1 = -f(r c )dt 2 + r 2 c dn 2 
= -dx 2 Q + r 2 c dQ 2 

= -^dr 2 + r c W, (7) 

where r = Xx° = Xy/Jt [14]. The aim for introducing the rescaling parameter A is 
to enable us to take the non-relativistic limit with ease: the limit just corresponds to 
taking A — > 0. Since the above line element contains no dr terms, we can think of the 
matrix representation of the corresponding metric tensor h^ v as containing a raw and 
a column full of zero corresponding to the r direction. Using (6), it is easy to check 
that h^y is idempotent, so this object can also be viewed as the projection operator 
onto the hypersurface. 

The extrinsic curvature (the second fundamental form) of the hypersurface is defined 

as 

^/M/ ~^£"n^ J ^iv (8) 

The projection of the Einstein equation onto the hypersurface yields some constraints 
over the geometric quantities on the hypersurface, and these will be best expressed 
using the first and second fundamental forms, as will be worked out below. 

There are two different types of projections, i.e. the longitudinal and normal pro- 
jections. These are given respectively as follows: 

(R^ - ^Rg^)n»h\ = (-Ag^ + T^)n»h v bl (9) 
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{R^ - ^Rg^)n»n v = {-Ag^ + T^nW . (10) 

These two equations are usually referred to as momentum and Hamiltonian constraints 
respectively. 

It will be convenient to make use of the Guass-Coddazi equation to turn the above 
constraint equations into a neat form. For this we recall that the following relations 
hold, 

A%, = hP a h p a h\h s v R a M& + {K\K m - K" v K aii ), 
DfrK^ = h\R pa n p , 

where R p a pv is the curvature tensor on the hypersurface and is the covariant deriva- 
tive acting along the hypersurface, D^h up = 0. Using the relation h ti l/ n 1 ' = and the 
Guass-Coddazi equation, we can calculate in a few lines to get the constraint equations 
on the hypersurface. The result reads 

D a {K\ - h\K) = T^nWt, 

R + K ab K ab - K 2 = 2 A - 2T^n fi n v . 

We shall see that these constraints reduce to Navier-Stokes equations on the hypersur- 
face if appropriate boundary conditions are imposed. 



3 Boundary conditions on the hypersurface 

The boundary conditions we would like to impose on the hypersurface is the so-called 
Petrov-like condition [20]. To describe such conditions, we need to project the bulk 
Weyl tensor onto the hypersurface. Inserting the Einstein equation (1) into the (n + 2)- 
dimensional Weyl tensor 

2 2 

Cni/fjp Rpuap {g 'a\aRp]i/ 9i/\aRp]p) H / ' 7T R 9 'u[cr 9 'p]v j 

n n[n + 1) 

we get the following relation, 

Cpvap = Rpucrp / : TT9p:[a9p]v {.9tAcrTp\ v gv[aTp\p)^ (H) 

n(n + 1) n 

where T pv and T are respectively the bulk stress-energy tensor and its scalar trace. 
Later on, we shall meet the following projections of the Weyl tensor, 

ri y.v y.<* up r* 

^abcd '» a"" 6'* c'* d^pvapi 

^abc(n) a'" fe" - c"- ^pvopi 

C a ( n )b(n) = h fJ, a n u h' T bn p C flU(7 p, 
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which can all be re-expressed in terms of the intrinsic and extrinsic curvatures of the 
hypersurface and the bulk stress-energy tensor, i.e. 

4(A - T) 

Cabcd = Rabcd ~ K ac K b( i + K a dK bc } — rrh a [ c hd]b 

n[n +1) 

2 

- -h p a h u b h a c h p d (g p[a T p]u - g^T^) , (12) 
2 

C a bc(n) = D a K bc — D b K ac — —h Pj a h u b h a c n p (g p [ <7 T p ] u — g v ^T p ^ p ), (13) 

Caln)b(n) = — Rab + KK ab — K ac K C b + a h° b Ru a ~, ' TT" h ab 

n{n + 1) 

2 

h p a n u h (T b n p (g p[a T p]v - g u[a T p]p ). (14) 

n 

The boundary conditions to be imposed on the hypersurface S c are given by [20] 

C{i)m = l^rnSFm/C^ = 0, (15) 

where the Newman-Penrose-like basis vector fields are introduced which obey 

I 2 = k 2 = 0, (k, I) = 1, {I, mi) = (k, rrii) = 0, (m*, rrij) = 6 l j. (16) 

In (n + 2)-dimensional spacetime, we have n basis vectors m„ i — 1, • • ■ ,n. is 
the /x-th component of m,. For our specific choice of spacetime metric (3), the basis 
vectors are given as follows [22], 

m = -— 1 —— — di, 

sin V\ ■ • • sin t/j_i 



J- (±r N 1 



fc = -Vliv7 S, + "J = "Vl (So + " ) ' 

Using these we can expand (15) as 

= 0, (17) 

where, by use of (12)-(14), we have 



CW)j — Roioj — K 00 Kij + K jK 0i — — — -/io[o^j]j 



4(A - T) 
n(n + 1) 
2 

-h p h 1/ i h cr h p j (g p[a T p]v - g u[a T p]p ), (It 



n 

2 

Coij(n) = D Kij — DiK 0j h p h u i h a jn p (g fJ ,[ rr T p ] l/ — g v [a-T p ]n), (19) 



n 

Ci(n)j(n) -^ij "I - KK^j K^K j -\- h p ih jR^u 12 (n I 1) 



- ^Xh^n^g^T^ - g v[tT T p]ll ). (20) 

In the absence of matter source, such as the case in [14], the dynamical degrees 
of freedom on the hypersurface are fully encoded in the extrinsic curvature K ab , or 
equivalently in the Brown- York tensor 

t ab = h ab K - K ab . (21) 

The total number of independent components of K ab is ( n + 1 K n + 2 ) _ However, the Petrov- 
like conditions provides ra( - n 2 +1 ^ — 1 constraints over the extrinsic curvature tensor, where 
the extra —1 is provided by the traceless condition for the Weyl tensor. Therefore, there 
are only n + 2 remaining degrees of freedom. These remaining degrees of freedom have 
to obey the Hamiltonian and momentum constraints described earlier in Section 2, 
which will be turned into Navier-Stokes equation for the fluid on the hypersurface. In 
this process, these n + 2 degrees of freedom will be turned into the density, pressure 
and velocity components of the boundary fluid. 

From the definitional) we get 

K\ = -h\t - t a b , t = t\ = nK\ = nK. 
n 

So, 

K = — K T = -t T ■ 

n 

n n 

Inserting these relations as well as (12)-(14) into (17), the boundary condition becomes 

2 t 2 t f t \ 2 k 

X 2 n 2 n \n J X y J ' 

- - {T vp n v n p - 2 A + T + Too - 2T p0 n p ) % + T tj - . (22) 

The calculation in getting this result is somewhat lengthy. We have made the coordi- 
nate transformation x° —> r, resulting in elicit coordinate index r appearing here and 
there, e.g. for the second and third terms on the right-hand side of (12), we have 

— K oo Kij + KojKoi 

= K\K l3 + K°jK i = K\K l0 + ^K\K T 3 

1 t 2 t t 

= f T f T ■ -I h - t T h - +t T t - t - 

X 2 n 2 n n 

where raising of indices is made using the metric h pu . We also used Roioj = 0, which 
is implied by our metric ansatz. 
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So far, we haven't paid a single word on the boundary conditions to be imposed 
on the matter source. In the presence of matter source in the bulk, there are some 
extra degrees of freedom for the matter. We also need to impose appropriate boundary 
conditions on these extra degrees of freedom. This will not remove all degrees of 
freedom for the matter source. The remaining degrees of freedom would then play 
some role in the fluid equation on the boundary. Naturally, the role played by these 
extra degrees of freedom will be some external impact force [18] [19] or fluid source. 
In this paper, we will restrict our discussion to spherically symmetric spacetimes with 
a perfect fluid source. In this case, the most appropriate boundary condition to be 
imposed on the matter source will be fixing p(r), p(r) and u T to be constant on the 
hypersurface at r = r c . We will see that such boundary conditions do work well in 
establishing the gravity-fluid duality. 



4 Near-horizon limit 

4.1 Perturbative expansion of the Brown- York tensor 

In this section we shall make perturbative expansions of the Brown- York tensor and 
estimate the order of various quantities that appear in the boundary conditions and 
the Hamiltonian and momentum constraints. 

The perturbed Brown- York tensor can be written as [14] 

oo 

t\ = Y,^ a t\ (23) 

n=0 

where we have intentionally taken the expansion parameter to be identical to the scaling 
parameter A appeared in (7), so that the perturbative limit A — > is simultaneously 
the non-relativistic limit. In this expansion, the n = term corresponds to the original 
(unperturbed) background. Using the definition (8), it is straightforward to calculate 
the components of the extrinsic curvature, which read 

K\ = ^7jd r f, K\ = 0, 

r 2yf f r 

Therefore, to the first order in A, we can write the Brown- York stress tensor as follows, 

r 
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where t = t a a is the trace of the Brown- York tensor. We assume that the smooth 
function f(r) has some zeros at r = r^, which corresponds to the event horizon of a 
black hole. In this case, we can expand f(r) near the horizon, i.e. 

f(r c ) = f'(r h )(r c - r h ) + ±f"(r h )(r c - r h f + • • • , 

where r c is the fixed constant determining the hypersurface S c . For the above expansion 
to make sense, we require the hypersurface to be located near the horizon, so that r c — 
is small. In particular, we let this difference to be related to A via r c — = a 2 A 2 , 
where a is an arbitrary fixed constant which is introduced to balance the dimensionality. 
Then, it is easy to see that provided 7^ 0, 



/ -> 0(X 2 ) + 0(A 4 ) + • 
- ^ £>(A 2 ) + £>(A 4 ) + 



V / 7= V(rc-r h )^/f'(r h ) + ±f»(r h )(r c 



r h ) + > C(A) + 0(A 3 ) + 



y^L ^ o(\) + o(\ 3 ) + ■ ■ ■ , 

r 

d r f = f'(r h ) + f"(r h )(r c - r h ) + ■ ■ ■ O(A ) + 0(A 2 ) + • • • , 

^^C( A -i) + 0(A) + ... . (24) 



4.2 Matter perturbations 



In the bulk, we take the matter source to be an ideal fluid, 

T/ja, = (p + P)un u " + P9i*v> ( 25 ) 

where p = p(r), p = p{r) because of spherical symmetry, and = (— -^,0, •■• , 0) 
is the proper velocity satisfying u 2 = —1. Recall that the bulk matter contributes 
in the boundary condition via (22). We can evaluate explicitly such contribution by 
calculating the relevant terms in (22): 

- -(T vp n v n p - 2 A + T + T 00 - 2T p0 n p )h ij + T tj 

= -l(-2A + 2p+(f-l)p)h ij . (26) 

Since we have taken p,p and u T to be fixed on the hypersurface, the only possible 
perturbation on the bulk matter must be realized via perturbing the other components 
of the proper velocity, i.e. 

u r = + Xu r{1) + ■■■ , 



u { = + Au l(1) + • • • . (27) 

Using = g^u", we can get 

u,= (y]+\u r ^ + --- , ^(A^W + •••)) • (28) 

In the next section, these perturbative expansions will be inserted into the perturbed 
version of the Hamiltonian and momentum constraints to yield the correct fluid equa- 
tions on the hypersurface. 



5 Fluid on the hypersurface 

Now we are in a position to combine all previous analysis to write down the fluid 
equations on the hypersurface. To do so, we first consider the perturbation of the 
boundary condition (22). With the aid of the order estimation (24), we get in the first 
nontrivial order O(X^) the following relation, 

a 3 3 u k) na 3 

- I (-2A + 2p - p + ^| r ,) hf - R v (29) 

where f' h represent the derivative of function / evaluate at and h[j is the metric of 
a standard n- sphere of radius r^, hf^ = r\dTL 2 n . 

Now, consider the momentum constraints on hypersurface. Inserting the definition 
of Brown- York stress tensor, we get 

D a t\ = -7>". (30) 

Let us expand this equation in components. First look at the r component. The 
left-hand side of this component reads 



D a t\ = D T t\ + Drf 
1 



r 



= D T t T T — —Di(t T jh* 3 ). (31) 

We can easily make some order estimation for each term appearing on the right-hand 
side, 

D T t\ = D T t T ^ ] + >0(A (1) ), 

^Diifjh 1 ' 3 ) = hi^Dff + -..-► O(X^). (32) 
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Therefore, we determine that the leading order of (31) is jh^^ D i t T ^ 1 \ The right-hand 
side of the r component of (30) is 



T^n" = T TT n T + T rr n r = (p + p)(«/))i T -»■ £>(A (1) ) + ■ ■ • . (33) 



Since the leading order of (p + p) is O(X^), u^n^ = -jjU^ + ■ ■ ■ is O(X^), and 



u T = a/7+ Aw r W + ■ • • is O(X^), the entire term is O(X^). So, we get the first useful 
equation at (9(A ( -~ 1 - ) ): 

Dit Tl{1) = 0. (34) 

Next we consider the spacial components of the momentum constraints. We have 

D a t a i = D T t\ + Djt j i. (35) 

Inserting (29) into the above equation, we get in the order O(X^) the following ex- 
pression, 

D T ff» + ( 2rf ) J D fc ri 1) + tfMDjfV - D k {Dit T W + D k ef ] ) + ^D k t^ - D,W r 

(36) 

It can be shown that the components of the Christoffel symbol on the hypersurface 
S c are vanishing, we have 

DrfV = d T t T f ] - f a Ti t T ^ = d T t T f\ (37) 

This gives a simplification to the first term in (36). The second term in (36) is zero 
thanks to (34). Using the same equation (34), we can also get 

D k D t t Tk ^ = (D k Di - DiD k )t Tk V = Ri m t rm{l) . (38) 

This quantity has appeared in the 4-th term in (36). Now, according to Bianchi identity, 
we have 

D a R a i = —DiR. 

Our specific choice of line element on the hypersurface implies that R T i is vanishing 
and R is constant. So, we have 

DaR 1 , = Dj&i = ^D,R = 0. (39) 

So, the last term in (36) is also vanishing. In the end, (35) becomes 

D a t\ = drfV + -^= (wMDjfV - D k D k t^ - R m l t T gA + (40) 
V h > n 

in the leading order O(X^). 
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Next, let us consider the contribution of bulk matter to the momentum constraint. 
Evidently, 



T Mi n M = T Ti n T + T ri n T = (p + p^u^n^Ui. 
The leading order of this term is O(X^), 

2>" = A(p + p)(tiX) (0) ^ 0) ii« i(1) + °( X(2) )- ( 41 ) 
It is crucial to introduce the following new notations, 

t T \ 1] = ~, — = K h = -T, l n^\ (42) 

2 n 2 M v 7 

Here ^ and p are to be interpreted respectively as the velocity and pressure of the 

fluid on the hypersurface. Inserting these new notations into (40) and (41) and then 
substituting the results into (30), we get the following equation, 

drV . + Dkp + " (vj D . v . _ DkD k v . _ Rm. v Jj = U (43) 



Recall that a is an arbitrary finite constant which is introduced to balance the di- 
mensionality in the relation r c — = a 2 \ 2 . We can freely choose a = \pf T h as did 
in [17]. Doing so (43) becomes the standard forced Navier-Stokes equation in curved 
background, i.e. 

d T Vi + D k p + v j DjVi - D k D k v,i = fi + R m iV m . 

We intentionally moved the curvature term to the right-hand side, because both fi and 
the curvature term provide external forces to the fluid equation. The force fi is related 
to the bulk stress-energy tensor explicitly via (42), while the curvature force is also 
related to the bulk stress-energy tensor implicitly through the bulk Einstein equation. 
The equation (34) is re-interpreted as the incompressible condition 

D k v k = (44) 

of the boundary fluid. 

While completing the draft of this paper, a new paper on gravity-fluid duality 
appeared on arXiv [23], which contains some overlap in interests with our work. The 
difference lies in that, in [23], gravity-fluid duality for general non-rotating black holes 
are established, while in our work only spherically symmetric black holes are considered. 
However, our work still stands in that we have included an ideal fluid source in the 
bulk, while in [23], only sourceless gravity is considered. Both works aim to extend 
gravity-fluid duality to more general setting. One can expect that more fruitful result 
in this field is upcoming. 
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